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!♦ Introduction 



Transonic flow in turbomachinery has become a field of increasing interest 
since requirements for higher engine thrus t-to-weight ratios have forced the 
operating conditions of compressors and turbines into the regime of mixed 
subsonic - supersonic flow conditions. Since transonic flows are usually 
accompanied by shock waves which may adversely effect compressor/ turbine 
efficiency due to wave drag, viscous interactions and unsteady effects, it is 
important for the design aerodynamicist to have computational tools for tran- 
sonic airfoil-, cascade-, and channel flows which minimize the shock effects. 
This is analogous to the situation in aircraft design where tools for wing and 
body definition to obtain shapes with minimal drag for given lift are needed. 

The flow in turbomachinery is usually unsteady because of blade row inter- 
actions. In addition, viscous effects are often present. Current analyses 
usually neglect these effects in order to arrive at models which are tractable 
by current computational methods. Hence, if we restrict our model to inviscid, 
steady, irrotational compressible flow past a blade row as shown in Figure 1, 
then we have at least a clear formulation of a defined goal: we intend to 

give shapes for the blades which will lead to shock-free flow with the above 
described idealizations and for special operating conditions. Figure 1 is an 
’’artist’s concept” of a flow with local supersonic regions and isentropic re- 
compression; sonic surfaces with a local ilach number M equal to One are 
drawn, dividing the supersonic regions (M > 1) from the surrounding subsonic 
field (M < 1) . A further idealization is the investigation of plane flow 
examples instead of the sketched three-dimensional problem. The unwrapping 
of cylindrical sections of a rotor or stator leads to infinite sequences of 
the airfoil-like blade sections in a two-dimensional plane. This - in compari- 
son to real flow in turbomachineiy" - modest configuration is our goal here. 
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The following chapter will introduce the usual basic equations for this 
type of problem. 

Following this a hodograph transformation and a method of characteristics 
will be presented. This chapter serves to explain a new mathematical princi- 
ple for mixed flow computation, which was first applied to solve hodograph 
boundary value problems. 

In the next chapter a direct method in the flow plane is developed using 
a more physically defined model rather than a mathematical principle. Both 
approaches are, in fact, identical procedures, but the latter direct method 
certainly will have more practical value, since a high amount of experience is 
usually needed to operate hodograph methods in practice. Nevertheless, impor- 
tant insight into and understanding of transonic flow phenomena may be gained 
from indirect methods and problem formulations. 

The presentation of design methods is followed by some examples obtained 
with these techniques. The direct design method is already operational for 
aircraft wing sections, therefore an example of direct shock-free single air- 
foil design is included. Development of single airfoil design procedures but 
also an extension to three-dimensional flows is given in Ref. 1. 

2. Basic Equations 

We consider steady compressible isentropic flow with the basic relations 
of continuity and irrotationality : 

div (pw) = 0 (la) 

rot (w) = 0 (lb) 

with the flow vector w and the density p . The wellknown gasdynamic rela- 
tions for isentropic flow define p , w , the speed of sound a and the Mach 
number H : 



2 



W = I w 



2 

a 



2 



a 

o 




M = 'W/SL 

p/p _ = (1 + M^) 

o Z 



1/Y-l 

= F(M) 



( 2 ) 



with Y the ratio of specific heats and subscript o defining stagnation 
conditions. For two-dimensional flow in an x, y-plane the relation (la) 
permits introduction of a streamfunction ; irrotationality (lb) allows 
definition of a velocity potential , with its gradients equal to the velocity 
components in x- and y- direction: 

i ^y " ^ ^ cos6 

(3) 



The flow angle is 6 . The system (3) for (f> and 'p is a generalization 
of the Cauchy-Riemann equations, elimination of ^p by differentiation yields 
a Poisson equation for 6 : 



P P 

<p 4* cb = - (-2S 0 + ^ 

^xx yy p p ‘y 



(4) 



with p a function (2) of M , and therefore 

p = p(w) = P(4>^^ + , (5) 

the system is nonlinear. Furthermore, the system is of elliptic type when 
M < 1 , and of hyperbolic type when M > 1 , with a parabolic type line along 
the sonic velocity locus M = 1 . 
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3> An Indirect Method to Compute Shock-Free Transonic Flows 



The following transformation into a special kind of hodograph variable is 
included to gain advantages of the resulting system against equation (4) : 
first, the equations will become linear, but also the sonic locus M = 1 will 
be known in advance so that the type change from an elliptic to a hyperbolic 
solver routine will not have to be iterated. We introduce a Mach number - 
dependent variable v(M) , which is defined differently for subsonic and 
supersonic Mach numbers; ** 



V(H) 



C arctan (C 



- 



1 ) - 



y 



C artanh (C v/l - M^) - artanh 



M^ - 


1 , 


for 


M > 1 


YU 


Y , 


for 


M £ 1 



( 6 ) 






xn 



Another variable used in the following relations is 

o\ 1/Y-l 



K(M) = K(M(\))) = ^1 






M - 1 



With the flow angle 



6 = arctan /u , 



(7) 



( 8 ) 



a working plane (v,6) , the socalled Rheograph - plane, is defined. As 
illustrated in Figure 2, the transformed system (3) in the Rheograph plane is 
linear but different in the supersonic (v ^ 0 , M ^ 1) and subsonic 
(v ;< 0 , M £ 1) half-plane, '^The coefficient K(v) is positive, with a 
weakly singular zero at the sonic condition V = 0 , 

The system for subsonic flow gives a linear (elliptic) Poisson equation 
for (p . 

%v ^ 10) , (9) 
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and a linear hyperbolic system for supersonic velocities: 



'^vv *^65 



(V > 0) ; 



( 10 ) 



The sonic locus v = 0 is included in both systems as a contact line 
with common sonic values 



<^^( 6 ) = (J) (V ^ + 0 , 6) , 
= Ip (v + 0 , 5) . 



( 11 ) 



In order to obtain any solution with these common values we will have to 
solve the two systems in some suitable coupled way. This seems difficult, 
but the following simple idea allows a decoupled computation, which will prove 
effective for design computations. 

A flow example with a local supersonic region embedded in subsonic flow 
can be obtained in the following way. At first we will solve a boundary value 
problem for the elliptic system (v ^ 0) in order to get the subsonic part of 
the flow and to obtain proper sonic line values along the 5- axis. The latter 
is effectively carried out by an analytic continuation of the elliptic boundary 
value problem into the hyperbolic half plane, Figure 3a , and solution of the 
elliptic system within these boundaries; the hyperbolic system (10) for v > 0 
is made elliptic like (9) in the continuation region, temporarily by a simple 
change of one sign. 

Assume we found a computational solution for this problem by use of an 
elliptic solver routine. The solution will be valid for local subsonic flow, 
values $*(6) =9(0 , 6) are avila le. A local expansion of the solution 
near V (-0) by use of the system in figure 2 and the structure of K(v) 
near v 0 , namely 
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( 12 ) 



K* (V) = ^ 





has to have the structure 




3 ^ * + 0 (v^) (13) 

6 



With a known solution (p in v ^ 0 the streamf unction ip*(6) along the sonic 

% 

line can be integrated now. 

Thus both potential and streamfunction are known along the line M = 1 , 
they form initial conditions for the following integration of the hyperbolic system 
in the halfplane V > 0 , (Figure 2) . 

Integration of the hyperbolic system is carried out by a method of charac- 
teristics which will be outlined here because of its use also for the direct 
method described later. 

Introduction of characteristic variables 



which means, that the inclination of characteristics C = const,, r) = const* is 
defined in the (j) ,i/;- plane: 




(14) 



(see Figure 3b) leads to a transformed hyperbolic system 




(15) 





( 16 ) 
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With known values ({)* , ij;* along the initial (sonic) . line, this relation 
(16) allows a marching procedure away from the initial l±n^. Figure 4 illustrates 
how the solution in an infinitesimal triangle is obtained: and ^ in 

point 3 is defined by given values in points 1 and 2 and the averaged coeffi-' 
cient K 

m 



Given: , 62 , , '1^2 ’ ^2 ’ ’^2 ' 

Asked: > ' 1^3 • 



^3 2 ~ *^ 2 ^ ’ 









(17) 



4) + (j) 4 ^ - 4» 

^3 = ^ • 



^ 1 + 4^0 , <{>■,+ 

ih = + -= = 

^3 2 K 2 

m 



Due to the analytic structure (12) of K near the sonic condition, for 
sonic initial points 1,2 the averaging of K is 



K (V. = V = 0) 
ml 2 



i «" 3 > 



(18) 



This ”sonic” averaging is a discretized description of the fact that the 
characteristics in the (}),i|;-plane, (as well as in the physical plane x,y) , 
have Neil’s parabola cusps when they meet the sonic line. Figure 5a . 

A mathematically possible, but physically unrealizable solution exists if 
the characteristics of one family (C or r| = const,) intersect each other, 
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Figure 5b, Though being useful for interpretation of shock-formation in a free 
flow, such solutions do not have much practical use. 

For the given initial conditions, 

V* = 0 , 6 > 6 > 6 : 

C{) = cti* (6) , (19) 

•» 

I}; = (6) 

the whole triangle ABC, Figure 4a, can be evaluated by a straightforward integra- 
tion based on (17) using the characteristic grid and marching into the V-direc- 
tion. This way, with the starting sonic isotach V- = 0 , neighboring isotachs 

+ Av = + Y A6 (20) 

can be obtained by this method, a concept which was extended by the author also 
for the computation of three-dimensional flow fields, Px.ef. 2. The solution 
(f) , Jp (v, 5) within the triangle ABC defines now also the solution in the 
physical plane x,y by 

’ 

( 21 ) 

“ ~7~~\ (sind d(J) + cos5 — ) , 

^ w(v) ^ p(v)^ 

particularly the shape of the resulting streamline il;(v,5) = , see Figure 2b, 

if the elliptic boundary value problem solved before (Figure 3a) was defined by 
a curve where jp - jp^ - const , which is a Neumann type boundary value problem 
for system (9). In Figure 3c the solutions of the elliptic (v _< 0) and the 
hyperbolic (v 0) region are combined now and, though solved by two different 
systems, they represent one transonic flow solution, which can be evaluated in 
the physical plane x,y by (21). 



8 



As we know from operational hodograph methods, Refs. 3 and 4, the rela- 
tions between the hodograph (input) boundaries and the physical flow (output) 
boundaries resulting from the obtained solution are usually complicated and it 
requires rather much experience to operate such indirect methods in order to 
obtain useful examples of transonic airfoils, cascades or other flow boundaries. 
Some ideas about airfoil design with methods of this type are given in Refs. 

5 and 6, here we only want to point out the basic idea of solving an extended 
elliptic boundary value problem first and use the result to integrate the 
hyperbolic stystem in the second step. This technique will be used again in 
the method outlined in the follow’ing chapter. It will be explained without 
using the details of the hodograph method, but the reader undoubtedly will 
find the same basic idea resulting in strong relations between the indirect 
and the new direct method. 

4. A New Direct Design Method for Transonic Configurations 

The isentropic relation (2) for the density 

p/p = F(M(({) ^ + * “)) (22) 

o x y 

is responsible for both the nonlinearity of (1) and the elliptic/hyperbolic 
mixed type of the Poisson equation (4). 

The following idea will not remove the nonlinearity of the system but it 
will establish a uniform elliptic type Poisson equation throughout the whole 
flow field and Mach number range. 

Imagine a modification of the isentropic gas model in the way that the gas 
might behave like an isentropic compressible gas as long as the local Mach 
number is subsonic. When the local speed of sound is exceeded in the flow 
past a given configuration, the density of the gas model is frozen to the 
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value of sonic flow condition. This means that this fictitious gas is incom- 
pressible in supersonic flow conditions! More generally we might think also of 
other relations p(M) when M _> 1 : 

p/o = F(M) : 
o 



F(M) = 



/F(M) 

\f(m) 



(M < 1) , 

(M > 1) , 



(23) 



F(l) = F(l) 



but we have to ask for possible constraints for F(M) in order to get an 
elliptic and not again a hyperbolic system similar to the original isentropic 
one. This constraint exists and it is defined by the limiting conditions 



p * w 
p • a* 




(M) 



M 



Y + 1 



a/2 



a + (y-i)m 



= const. 



yielding 



with 



F (M) = F’^ 



2 4 (y-I)M^ 
(Y+1)M^ , 



1/2 



F* = p*/p 

o 




(24) 



(25) 



(26) 



A fictitious gas with F = F^ in the supersonic region will lead to para- 
bolic equations within this region. This model is therefore a limiting and 
excluded border of a whole area of possible relations for F , see Figure 6 . 

A one-parametric family of possible fictitious elliptic gases is defined by 
the density - Mach number relation for ^ 1 
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Fp(M) = F* 



' 2 -I- (Y-l)M^ 
\ (y+1)M^ 



P/2 

I . 



( 27 ) 



P < 1 



excluding the parabolic limit P = 1 and including the above mentioned "incom- 
pressible supersonic flow" when P = 0 , see Figure 6. 

The isentropic relation F(M) and the parabolic limit have the 

same gradient at M = 1 , thus indicating that the isentropic relation leads to 
parabolic equations just at the sonic surface in the flow field. 

With a p(M) relation (23) the basic Poisson equation (4) will be non- 
linear but of elliptic type throughout the flow field: 



(J) + (j) = P 

XX yy 




(28) 



Corresponding to the real and the fictitious part of the density relation 
(23), a solution ^(x,y) will consist of an isentropic part cf) and a fictitious 
part i , see Figure 7a : 



^4(x,y) 

~ / 

4> = 



where 



where 



<b ~ ^ ^ 

X y — 1+Y 

~ 2 -.2 2 a 2 



(29) 



but with common values <J) and (p along the sonic line 



(J)(x*,y*) = ^ (x*,y*) = 4)* , 



where 





(30) 



From a continuous fictitious density relation (23), namely F(l) = F(l) , 
follows that not only (b but also grad ^ is continuous across the sonic line. 
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Thus the velocity components along the sonic line may be given by 



(u* , V*) = grad (}) = grad (}) (x* , y*) = (<}) 



X 






( 31 ) 



If we define now all sonic condition data as functions of the sonic line arc 
length s , with 

2 2 2 
ds = dx* + dy* , 



X* = x*(s) , 

y* = y*(s) , 

u* = u*(s) 

V* = v*(s) , 

6* = 6'*(s) = arctan '' /u* , 

(p* = <p*(s) = Cfi (x*(s) , y*(s)) , 



(32) 



V* , 



then the s treamf unc tion 'p* along the sonic line follows from (21) to 



(|)*(s) = 



r • 

• F* • J (u* y* “ V* X 



X*) ds . 



(33) 






With no sources within the local fictitious supersonic region, which is 
bounded by a wall and the sonic line, see Figure 7a, the integral over the 
sonic line arc AB has to vanish: 



ip*(0) = = 0 • (34) 

This can be used, if needed, to check the accuracy of the conservative 
numerical algorithm for the elliptic system solved before* Numerical inaccuracy 
gives a result , in which case a refinement of the elliptic computa- 

tion or some correction should be made in order to obtain conservation of mass 
implied by (34) up to a sufficient degree of accuracy. This having been 
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achieved, we have from (32), (33) the variables 6* , (j)* , 4;* , x* , y* as 
functions of the parameter s which we will now need as initial conditions to 
solve the real isentropic supersonic flow domain and redefine the wall stream- 
line = 0 between the points A and B in the physical plane. Since (32), 

(33) is equivalent to (19), the method of characteristics as described above 
for the indirect approach can be used again. Figure 7b shows a sketch of the 
result, the characteristics grid and the new surface shape which is compatible 
with the computed shock-free transonic flow. Compared to the original shape, 
the modified one is flatter providing more space for the real supersonic flow 
to pass than the fictitious flow with higher density (F > F I) required before. 

It can be shown that the streamlines of the obtained hyperbolic flow 
solution fit into the subsonic flow x^ith smooth slope and curvature. This is, 
of course, important for the practical value of the obtained new wall shape, 
see Figure 7c. The possibility of obtaining limit lines was illustrated in the 
previous chapter. Figure 5b; the mentioned restrictions for practical use of 
such results apply of course also to the approach outlined here, since the 
method of characteristics is used in the indirect as well as in the new direct 
me thod. 

The following chapter will present some examples of shock-free transonic 
flow through cascades and past airfoils obtained with the indirect (hodograph) 
and the direct design technique which were both discribed here. 

5. Application to the Design of Transonic Cascades 

The indirect method described in chapter 3 was at first carried out by 
means of an analog elliptic solver. The rheoelectric analogy was used to solve 
the elliptic boundary value problem in the Rheograph (v,5-) working plane. 

This analogy has the advantage of a physically clear and illustrative formulation 
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of the type of problems studied, A systematic technique of obtaining super- 
critical airfoils was developed, see Ref. 6, 

Despite interesting results which led to the development of practical 
supercritical wing sections, the method could hardly be put into operation on 
a broader scale because of the electromechanical work involved with the set-up 
of boundaries, measurement, etc. The analogy has been successful in many fields 
not only in fluid mechanics but the fast development of larger and faster 
digital computers opened other, and in many respects more preferable, ways to^ 
compute problems where the analogy has served before. For use on a wider scale 
in the aircraft industry the analog set-up was, therefore, replaced by computa- 
tional elliptic solver routines like panel methods. Ref. 7. The most recent 
development is the use of the fast Poisson solver for the elliptic part of the 
problem. 

A result of this more economical way by using numerical techniques instead 

of the analogy is shown in Figures 8-11. Figure 8a depicts the elliptic boundary 

value problem in the V, 6-plane , with singularities, a second Riemann^s sheet 

and a branch-point , but also with an elliptic continuation region reaching into 

the hyperbolic half-plane V > 0 . The complicated structure of the elliptic 

boundary will not be discussed here, it defines the resulting physical flow. 

In this case a compressor cascade was designed, as illustrated below. Figure 

8b and Figure 9 show the results of the hyperbolic computation using the method 

of characteristics described in a previous chapter. In Figure 9a the resulting 

characteristics grid in the (J),i|;-plane is drawn, in Figure 9b the same 

characteristics are shown in the physical plane with the streamline \p = 0 . 

o 

In Figure 9c the isotachs V(M) = const, and the isoclines C = const. are 
drawn, thus defining the whole supersonic flow field. The resulting cascade 
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is drawn in Figure IQ , the lower blade is the inviscid design with the local 
supersonic flow field. It has an open trailing edge to allow a boundary layer 
displacement subtraction. A boundary layer computation was made using the 
design pressure distribution depicted in Figure 11 . The resulting actual blade 
is shown in the upper part of Figure 10. The pressure on the rear part of the 
blade follows a Stratford distribution, which is favorable in avoiding boundary 
layer separation. The thickness of this blade has a minimum somewhat upstream 
of the trailing edge, which results from the chosen elliptic boundary in order 
to obtain this type of pressure distribution, see also Ref. 8. 

The example of this compressor cascade design was obtained by the indirect 
method. As outlined, the direct approach has the advantage that an airfoil- or 
cascade- geometry can be chosen and only the part wetted by supersonic flow 
will be modified by the algorithm. In Ref. 1 a new combined analysis/design 
method for 2D airfoils is described. It is based on Jameson’s circle plane 
computer code for analysis of airfoils in transonic flow. Ref. 9. In Figure 12 
a conventional NACA 64A410 section was modified by this procedure in order to 
obtain shock-free flow past a single airfoil. The pressure distribution on the 
NACA airfoil (with a recompression shock) and on the airfoil obtained by the 
design procedure (shock-free) are compared. Little lift change and vanishing 
drag are the practical results of this inviscid computation. Of course, viscous 
effects can be introduced here also by starting with an open trailing edge 
airfoil shape. 

In a magnified scale the airfoil geometries are compared illustrating the 
”shave-off” on the upper side. Similar results for compressor cascades like 
the one in Figure 10 could be obtained also by the direct method if the input 
blade were thicker in the region which came out supersonic by the indirect 
method. This is sketched in Figure 13 . The implementation of the direct 
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design method into a cascade computer program is currently under way and will 
be reported in a forthcoming paper. 

A final remark concerning the design of accelerated flow turbine blades 
should be made here: the indirect method was used for design of the transonic 
inlet of a cascade with flow accelerated to supersonic exit velocity, see 
Figure 14 b,c . A result of this type can be obtained by the direct method if 
the chosen input cascade does not have closed blades. Figure 14a. Streamline 

lb 

contraction of the fictitious "elliptic supersonic flow" due to blade displace- 
ment will provide the necessary acceleration through sonic conditions. In 
Figure 14b the method of characteristics starting at the computed sonic line 
provides the transonic part of the field M > 1 , and finally a continuing 
method of characteristics using data on the limiting characteristic BC and 
given wall geometry BT", CT' closes the designed cascade blade. Figure 14c. 

Here the contour part AC of the blade is given by the algorithm. Of course, 
the streamfunction distribution 4 ^* may not have equal values , as 

postulated in (34), but 



^B “ ^A " 



case 



(35) 



case 



with Aij; 



being the amount of flow passing between two neighboring blades. 



6, Conclusions 



A new idea to compute transonic flows in turbomachinery has been outlined 
and demonstrated as the basis for computational codes to obtain shock-free flow 
boundaries like airfoils and cascades. An indirect method, which was developed 
first, served to outline the principle of elliptic continuation as well as a 
method of characteristics. Both techniques were then used in a direct approach 
to design shock-free configurations. This method can be used in combination 
with any flow analysis code thus extending such a program to a combined 
analysis/design tool. 

Conventional airfoil or cascade geometries can be made shock-free by 
relatively small alterations of the shape. In combination with a viscous dis- 
placement concept boundary layer effects may be included and will not be 
complicated for design operating conditions by shock - boundary layer interaction. 
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1. SHOCK-FREE TRANSONIC FLOW IN TURBOMACHINERY 
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2. POTENTIAL FLOW DIFFERENTIAL EQUATIONS IN THE RHEOGRAPH PLANE 
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ELLIPTIC CONTINUATION IN THE RHEOGRAPH PLANE 





4. METHOD OF CHARACTERISTICS, STARTING AT THE SONIC LINE 





5. CHARACTERISTICS NEAR (a) SONIC LINE AND (b) LIMIT LINE 
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6. ISENTROPIC AND FICTITIOUS RELATIONS BETWEEN DENSITY AND 
MACH NUMBER 
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7. ELLIPTIC CONTINUATION IN THE PHYSICAL PLANE: (a) 

Srwure resulting WAL 
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FICTITIOUS GAS 
SHAPE WITH 
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8. RHEOGRAPH BOUNDARY OF A SUPERCRITICAL CASCADE FLOW: (a) ELLIPTIC; 
(b) HYPERBOLIC PROBLEM 



26 



PSI • 

0 . 16 - 






9. COMPRESSOR CASCADE DESIGN: LOCAL SUPERSONIC FLOW FIELD (a) 

plane, (b) x,y-plane (CHARACTERISTICS), (c) x,y-plane 
(ISOTACHS AND ISOCLINES) 
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M^-0.S7 



10. COMPRESSOR CASCADE DESIGN: LNVISCID DESIGN, SUPERSONIC REGION, 

BOUNDARY LAYER SUBTRACTION (Re = 1 , 6 . 106 ) 
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n. COMPRESSOR CASCADE DESIGN: PRESSURE DISTRIBUTION 
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1.20 0.80 0.40 - 0.00 - 0.40 - 0*80 - 1.20 




NfiCfl 64fl410 mCH = 

ORIGINRU 0 0 o ) 

CL 3.6708 
CO 0.0032 • 

CM - . 1 489 



720 ALPHA =0.0 

OPTIMIZ£D( j 

0.6281 

0.0000 

-.1385 




COMPARISON OF THE PRESSURE COEFFICIENTS, SONIC LINES AND SHAPES 
FOR THE BASELINE NACA 64A410 AND THE SHOCK-FREE AIRFOIL OBTAINED 
FROM IT BY THE 2D DESIGN PROCEDURE 
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13. DIRECT DESIGN METHOD: ORIGINAL ( ), AND RESULTING SHOCK-FREE 

( ) CASCADE DESIGN 
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14. TURBINE CASCADE DESIGN: (a) ORIGINAL CASCADE, (b) TRANSONIC 

INLET, (c) RESULTING CASCADE 
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